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ABSTRACT 

We divide a sample of 302 type-1 AGNs into two subsamples based on 
the narrow line [O III] /R^nlr ratio, since we expect that there will be a 
stronger starburst (H II region) contribution to the narrow line emission for 
i?=log([0 lll]/~H./3 NLR ) < 0.5. For the two samples we find significant differ- 
ences in correlations between spectral properties of objects with R < 0.5 and 
R > 0.5. We find similar differences when we divided the sample based on 
the FWHM ratios of [O III] and broad H/3 lines (i?i=log(FWHM[0 III]/FWHM 
H/3broad)< — 0.8), i.e. similar correlations between R > 0.5 and Ri < —0.8 sub- 
samples from one side and R < 0.5 and R\ > —0.8 subsamples from the other 
side. The most interesting difference is in the correlation between the broad H/3 
FWHM and luminosity in the R < 0.5 (R\ > —0.8) sample that indicates a 
connection between the BLR kinematics and photoionization source. We discuss 
possible effects which can cause such differences in spectral properties of two 
subsamples. 

Subject headings: galaxies: active - galaxies: emission lines 
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Introduction 



Spectral properties of Active Galactic Nuclei (AGNs) depend on physical conditions and 
geometr y of emission regions which are quite likely to be changing during AGN evolution 



(see e.g. 



Lipari fe Terlevich 



20061 : 



Wang fc Zhang 



20071 ). To understand the complex 



properties, trying 1 


;o fine 


and explain 


leir i 




lysical background 


(see e.g. ] 


3oroson & Green 


1992: 


Wills et al. 


1999: 


Croom et al. 


2002 




3hang et al. 


2003: 


Yip et al. 


2004; 


Grupe 


2004; 


Wang et al. 


2006 


2009: 


Ludwig et al. 


2009; 


Kovacevic et al. 


20K 


), 


and references 
properties in a 


therein) . 


Boroson & Green 


(11992 


) (BG92) found a number of AGN spectral 



sample of AGNs, as e.g. that as the optical to X-ray slope and equivalent width (EW) of 
Fe II increase, the EW of [O III] decreases, as well as FWHMjj of H/3. These correlations 
are part of the eigenvector 1 (EV1) in BG92 principal component analysis. The underlying 
physics driving EV1 is not yet completely understood. 



On the other hand, the Baldwin effect ([Baldwin 



19771 : 



Baldwin et al. 



19781 ) is present 



in a number of lines, i.e. their line EW decreases with increasin g continuum lumino sity. 



The majority of the UV and optical lines show a Baldwin effect (jDietrich et al. 



20021 ). but 



for some lines there is no trend with the continuum luminosity, or e ven an inverse Baldwin 



effect . The inverse Baldwin e ffect is found f or the Fe II optical 



2007 



2007m . 



Kovacevic et al. 



Kovacevic et al. 



2010h and for H/3 (lOroom et al 



2002 



i nes (INetzer fe T rakhte nbrot 



Netzer fe Trakhtenbrot 



(120101 ) found an inverse Baldwin effect for the broad H/3 component 



(in FWHM H/3b ro ad <3000 km s 1 subsample), but on the other hand, the narrow H/3 



component shows a normal Baldwin effect as do the [O III] lines. The origin o: 



f the 



Baldwin effect is an open question still under debate (for review see 



Green et al. 



2001 



and reference therein). Some indications of connection between the BG92 EV1 correlation 



x Full Width at Half- Maximum Intensity of the H/3 broad component 
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with the Baldwin effec 


t are discussed in 


the literature 
)■ 


Kovacevic et al. 


2010; 


Zhane et al. 


2011 



There are some indications that the EV1 may be connected with AGN evolution (see 



Marziani et al. 


2003; 


Grupe 


2004; 


Wane et al. 


2006 


), while evolution of AGNs is probably 


related with starburst regions Jlji'pari & Terlevich 


2006 




Mao et al. 


2009 




Sani et al. 





2010l ). Namely, it is possible that AGNs in an earlier phase of their evolution are composed 



of starburst (star-forming) regions and the central engine (A GN), and, during evo 



the starburst con tribution becomes weaker and/or negligible (IWang fc Wei 



Mao et al. 



20091 ). 



2006 



ution, 



2008 



Authors of recent papers (see e.g. 



Yip et al. 


2004; 


Grupe 


2004: 



Ludwig et al. 



2009J) 



used a large sample of SDSS (Sloan Digital Sky Survey) objects and applied Principal 
Component (PC) analysis on some spectral parameters. They noticed that PC projections 
depend on the sample properties; i.e., different divisions of sample by luminosity, [O III] 
strength, FWHM H/3 etc. lead to different correlations between spectral parameters, 
sometimes even opposite correlations for different subsamples. Some authors found that the 



significance of the Baldwin effect depends on the FWHM 



Sulentic et al. 



2009 



Zamfir et al. 



2010 



i0 ra nge of the sample (see e.g. 



Kovacevic et al. 



20101 ). 



In this paper we investigate a sample of the broad-line AGNs. Using the ratio of the 
[OIII] and narrow H/3 line we divide the sample in two groups and consider the connection 
between the BG92 EV1 and Baldwin effect into these groups separately. 

The paper is organized as following: in §2 we shortly describe the sample (subsamples) 
and methods of analysis, in §3 we present the obtained correlations and discuss the results, 
and in 54 we outline our conclusions. 
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2. The sample and analysis 



The sample o 



Kovacevic et al. 



302 AGNs used in this paper is the same as the one presented in 
( 120101 ). The sample has been chosen from the SDSS Release 7. It 



contains the broad- line AGNs within approximately uniform redshift range z=0-0.7, mainly 
within the lumi nosity range of 44<log(A Lsmn)<46. Other details of the sample selection 

J2010I). 



can be found in 



Kovacevic et al. 



As is described in 



Kovacevic et al. 



( 1201 Oh . we fitted the Fe II lines with a new 



template which enables precise estimation of the Fe II emission within the 4400-5500 A 
wavelength range (see Fig. [Q. The Balmer lines were fit with three components: a narrow, 
an intermediate and a very broad component (H/3 NLR, ILR and VBLR, respectively). The 
H/3 broad component is taken as the sum of the H/3 ILR and H/3 VBLR components (see 
Fig. E]). The FWHM and FWM10% (full width at 10% of maximum) of the broad H/3 line 
are measured as it shown in Fig. [21 

Both components of the [O III] AA4959, 5007 A doublet originate from the same lower 
energy level and both have the negligible optical depth since the transitions are strongly 
forbidden. We thus assumed that the [O III] A4959 A and [O III] A5007 A lines have the 
same emission-line profile; i.e. we fit each line of the doublet by one Gaussian (or by two, 
in the case of a significant asymmetry) assuming that the [O III] A4959 A and [O III] 
A5007 A lines have the s ame widths and shifts. We also fixed their intensity ratio at 2.99 



( IDimitrijevic et al 



20070 . 



To estimate the NLR contribution to the total H/3 flux (i.e., emission from the same 
region emitting the [OIII] lines) we introduce following assumptions: (i) only one narrow 
H/3 component is present in the total H/3 emission and (ii) this component has the same 
kinematical parameters as the [O III] lines, i.e. the same width and shift as the [O III] 
lines. Very often a blue asymmetry in the [O III] lines was present and in such case each of 
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the [O III] lines was fitted with two Gaussian components, but as a rule only one (central, 
non-shifted) has the same kinematics as the narrow H/3 line. 



To plot a "BPT diagram" ([Baldwin. Phillips fc Terlevich 



19811 ) shown in Fig. [31 we 



also estimated the narrow Ha component, connecting the width and shift of the narrow 
component with ones of the [Nil] lines. As it can be seen from Fig. [31 a number of AGNs 
from the sample show a significant starburst contribution. It is interesting that these AGNs 
are mostly with FWHM H/?broad < 3000 km s~ x (open circles), while in the AGN part of the 
BPT diagram both fractions (with FWHM H/3 broad < 3000 km s~] and F WHM H/3 broad > 



Sani et al. 



(120101 ) also found more 



3000 km s" 1 - full circles) are present. Note here that 
intense circumnuclear star-forming activity in Narrow Line Seyfert Is (NLS1) than in the 
broad-line AGNs. 

Since we have complete measurements of line parameters only for the narrow 
H/3 and [O III] lines in the whole sample (302 AGNs), we accepted a criteria of 
R=log([0 IH]/B.(3nlr) = 0.5 (horizontal dashed line in Fig. [3]) as an indicator of the 
predominant starburst emission contribution to the narrow emission lines. We divided our 
sample into two subsamples: R <0.5 (91 AGNs, hereafter "SB", starburst dominant) |^| and 
R >0.5 (210 AGNs, hereafter "AGN", AGN dominant). We should note here that such 
a high number of SBs is probably caused by the selection effect, since one of the criteria 
was that equivalent wi dths of typical absorption lines are small (for more details see in 



Kovacevic et al. 



2010). 



One could expect that where there is a signi ficant starburst co ntribution, the widths of 



the [O III] lines would depend on the luminosity (IBrungardt Ill988l ). We therefore plot the 



2 One spectrum (SDSS J130357.42+103313.50) is excluded from the sample because the 
[O III] emission is comparable to the level of noise. 
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log(L 5 ioo) against the log(FWHM [O III]) (see Fig. @| and find a good correlation between 
the continuum luminosity and FWHM of the [O III] lines in the SB sample R < 0.5 (r = 
0.70, P = 7.3E-15), while in the AGN sample there is no correlation. This indicates that 
even for the roughly estimated fluxes of the narrow H/3 component there is some difference 
between these two sub-samples. In Fig. [5] we plot the FWHM of the broad H/3 against the 
FWHM of the [O III] for all data, where the SB sub-sample is denoted with the full circles. 
As one can see, there is a good correlation between the widths for the SB sub-sample (r = 
0.74, P < 1E-16). 

It is a hard task to separate AGN s, HII regions and LINERs, even if one deals with on ly 



narrow emission line objects (see e.g. 



Rola et al. 



1997 



Dessauges-Zavadsky et al. 



2000 } . 



Therefore, the construction of BPT diagrams of broad line object s, where the permitt ed 



Nagao et al. 



(1200 lh 



lines are composed of two or more components, is very difficult, 
constructed the diagrams for a group of broad and narrow line Seyfert galaxies, and found 
that some Syl may be located in H II regions (see their Fig. 4), but there are no significant 
trends in the forbidden line ratios between NLSyl, Syl and Sy2 AGNs. Forbidden lines can 
be more precisely fitted than permitted lines, but even though forbidden line ratios depend 
on physical properties (electron density and temperature) it is hard to separate the objects 
using only forbidden narrow lines. 

Estimation of the narrow H/3 and Ha components (in the case of a broad line AGN) is a 
very complex task. There are several problems, and one is that in principle the line profiles 
(especially of the broad lines) are not pure Gaussians, i.e. , due to complex kinem atics in 



Smirnova et al. 



20071 ) a single 



the NLR (e.g., blue asymmetry produced by outflows, see 
Gaussian very often cannot represent the shape of the narrow H/3 and [O III] components. 
Consequently the approximation of two Gaussian functions could be used, which gives some 



indication about velocities of the subregions. Another problem in the estimation of the 
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narrow H/3 contribution could be that the solution is non unique, i.e., a combination of 
different Gaussians can give different contributions of the narrow Hf3 component. Therefore 
when fitting with multi-Gaussian components we are able only to estimate very roughly 
the contribution of the H/3 NLR component from the region emitting the [O III] lines. We 
estimated errors using a \ 2 test ( see Appendix I) and found that for the SB subsample, 
around 25% of the AGNs (those between estimated error-bars, R + AR >0.5) could belong 
to the AGN subgroup, but when we exclude these 22 objects, the correlations we obtained 
were not significantly changed. 

We also searched for another method to separate these two groups. We investigated 
the line width ratios of [O III]/H/?broad as a function of the continuum luminosity and found 
that there is no correlation (see Fig. [6]), but points with R < 0.5 (solid circles), mostly have 
i?i=log(FWHM [O III]/FWHM Hf3 broad ) >-0.8. Consequently, we also divided our sample 
into two sub-samples with Ri < —0.8 (hereafter the "AGN1" subsample) and R\ > —0.8 
(hereafter the "SB1" subsample). 

3. Results and discussion 

We calculated Spearman correlations between the spectral properties for the total 
sample of 302 AGN spectra, and for the "SB" and "AGN" sub-samples with R <0.5 and 
R >0.5 respectively (which correspond closely to R\ <-0.8 and R\ >-0.8). The results are 
given in Tables [D - @] and Figs. [7] - ED 

One of the most interesting correlations is between the FWHM (FWM10%) of the 
broad H/3 vs. co ntinuum luminosity (AL5100) (a similar correlation was found by 
Steiner I ( jl98lh for the FWZI of Ha) and the FWHM (FWM10%) H/3 vs. EW Fe II. 

These correlations are quite different for the SB and AGN subsamples. 
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The velocity-luminosity relation in AGNs has been discussed earlier. IShuder I (119841 ) 
found that t he FWHM of the Balme r lines increases with luminosity for a sample of 25 
AGNs. Also 



Wandel k Yahil 



(119851 ) reported about a correlation (r=0.5) between the H/3 



full width at zero intensit y (FWZI) and the 4 000 A continuum luminosity for a literature 



compilation of 94 AGNs. 



Stripe et al. 



(119991 ) found a weak correlation between the width 



of H/3 and lu minosity in the cas e of 1 26 AGNs. The same tre nd, i.e. a weak correlation 



was found by 



Joly et al. 



(I19851 ). but 



Boroson &: Green! ( 119921 ) found an anticorrelation 



(r=-0.275) between the H/3 FWHM and absolute V magnitude, significant at the 99% 
confidence level, in a sample of 87 PG quasars. They also noted that the FWZI of H/3 is 
quite sensitive to noise and depends strongly on the quality of the Fe II subtraction. We 
therefore measured FWM10% 

In Fig. 0we plot log(AL 5100 ) vs. log(FWHM H/3) and log(FWMI 10% H/3) of the 
broad H/3 component. We found a strong correlation between the continuum luminosity 
and FWHM H/3 for the SB and SB1 subsamples, while in the AGN and AGN1 subsamples 
(see Figs. [T] — [9]) , the correlation coefficients are much lower and not statistically significant. 
The correlation between continuum luminosity and FWMI 10% H/3 is present in the SB 
and SB1 subsamples, but not in the AGN and AGN1 subsamples (see Tables [1] and [2]) . 

The relationship between the EW Fe II and FWHM (FWM10%) H/3 is also different 
for the SBs and AGNs subsamples (see Fig. [8j Table [T|). Namely, the well known 
anti-correlation of EW Fe II and FWHM H/3, which is part of the BG92 EV1 correlations, 
is noticeable only for the AGN and AGN1 subsamples, while the SB and SB1 subsa mples 
show an opposite but statistically insignificant trend. Note here that 



Grupe 



(120041 ) also 



found a positive correlation between the EW Fe II and FWHM H/3 for NLS1 objects, and a 
negative correlation for broad line Syls. 



Some other correlations shown in Tabled] are also sample dependent: i.e., there is a 



- 10 - 



significant difference in the corrlation coefficients for the AGN and AGN1 than for the SB 
and SB1 subsamples. As can be seen, the Baldwin effect for [O III] and H/3 NLR tends to 
be stronger in AGN and AGN1 subsamples than for the SB and SB1 subsamples. 



3.0.1. Principal component analysis 

We applied principal-component (PC) analysis for the total sample, and subsamples. 
The results of the PC analysis for the two subsamples and for the total sample are shown in 
Tables [3] and HI As can be seen, the AGN dominant subsamples (AGN and AGN1) may be 
described with the first three eigenvectors which represent ~ 82% of variance (see Table [3] 
andHj). The first eigenvector, which describes ~ 39% of variance, represents dependence of 
line EWs vs. continuum luminosity (i.e., the Baldwin effect). The EW Fe II vs. EW [O III] 
anticorrelation is clearly present, implying that it is related to the Baldwin effect. Note, 
that in this eigenvector there is no influence from the FWHM H/3. The second eigenvector 
represents a relationship between line EWs and FWHMH/3, which should be similar to the 
BG92 EV1, but in our case, the EW Fe II vs. EW [O III] anticorrelation is not detected. 
In our second eigenvector we found the anticorrelation of FWHM H/3 vs. EW Fe II and 
the expected correlation between FWHM H/3 and EW H/3 broad. Our third eigenvector is 



dominated by the broad H/3 compo nent and EW Fe 
BLR and Fe II emitting region (see iKovacevic et al. 



I ind icating; a connection between the 



201 



The PC analysis results for the SB and SB1 subsamples are significantly different (Table 
|3]and|l]). THe variance may also be described with the first three eigenvectors representing 
84% (SB) and 76.5% (SB1) of variance. The first eigenvector accounts for 48.5% and 41.4% 
respectively of the variance of the subsamples, which implies that it dominates in both SB 
subsamples. It represents a relationship between the line EWs and continuum luminosity 
(the Baldwin effect), but also correlations/anticorrelation between the EWs with the broad 
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H/3 FWHM, since a strong correlation between FWHM H/3 and AL5100 is present. The EW 
[0 III] vs. EW Fe II anticorrelation is present also as the Baldwin effect of the lines. This 
eigenvector does not represent the BG92 EV1 since the FWHM H/3 and EW Fe II have a 
positive correlation. 



3.1. Possible physical interpretation of results 



The most interesting result is that SBs subgroups show correlations between the 
narrow and broad line widths and luminosities while the AGN subgroups do not. It has 
previously been noticed that some correlations between spectral propertie s depend on the 



FWHM H]3 range : 



2009 



Zamfir et af- 



file observed sample ( see 



2010 



Kovacevic et al. 



Steiner 



1981 



Grupe 



2004 



20101 . etc.). For example, 



Su 



Grupe 



entic et al. 



(120041 ] found 



that the FWHM H/3 and EW Fe II anti-correlate for AGNs with broad lines (BLSyls), but 
anti-correlate for NLSyls. Here we found two subsamples, one where there is a significant 
correlation between the FWHM H/3 broad component and continuum luminosity, and 
another where this correlation is not present. It is interesting that the correlation between 
continuum luminosity and FWM10% H/3 is the same as for the FWHM H/3 (see the SB and 
SB1 subsamples in Tables [T] and [2]) . 

Among other things, the different behavior of the two groups may be caused by: a) 
different contributions of starburst regions, i.e., a different nature of the ionization source, 
or b) by different Eddington ratios and/or geometrical structures of the BLR. 



3.1.1. A starburst contribution - a difference in the source of ionization? 



One can see from Tables [U and [2] that there are significant differences in spectral 
properties between the two subsamples. A possible explanation of these differences may be 
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that there is a difference in the source of ionization. 

As we mentioned above, there is a possibility that an AGN spectrum in the earlier 



(see 


Croom e 


; al. 


2002; 


Wane & Wei 


2006 


Popovic et al. 


( 


200 


9) found that in the case 



2008 



Mao et al. 



20091 ). Moreover, recently 



ratios correspond to starbursts rather than to an AGN origin. 

In the case of a dominant contribution from a starburst region to emission in the 
narrow lines, one can expect that, besides the AGN emission, additional significant starburst 



emission may be present in the composite continuum. Then, t here is probab 



source of continuum emission (as in the case of Mrk 493, see iPopovic et al. 



van extensive 



20091 ) that 



may be produced by a number of exploding stars. Consequently, the luminosity of such 
objects will depend on the star forming rate. In the starburst-dominant subsamples there 
is a positive, but statistically insignificant trend between the FWHM H/3 and EW Fe II, 
and it seems that in general the Fe II emission is stronger in this subsample (log(EW 
FeII)~1.7-2.3). The strong Fe II emission and weak [O III] emission ca n be explained with 



Lipari fe Terlevich 



a mo del that contains a massive starburst (SB) plus an AGN (see e.g. 
20061 ). The SB+AGN can lead to large scale expanding supergiant shells. A good correlation 
between the FWHM and luminosity in one group of AGNs may indicate that the broad 
component is not primarily a function of geometry (e.g., by rotation due to the gravitational 
force), but rather by random motion of the gas caused by different effects (such as the 
gas being randomly accelerated in several bursts). In such a case part of the flux of the 
broad Balmer emission lines may arise i n the stellar envelop es of Wolf-Rayet and OB stars 



associated with multiple SN events (see 



Izotov et al. 



2007 



and references therein). Note 



that EW H/3 broad only increases with AL5100 in the SB and SB1 subsamples, and this 
increase is not found in the AGN subsamples. 



13 



One can speculate that the FWHM of the broad component might not be formed as 
in the classical broad line region (where rotation velocity should b e present), bu t as a sum 



Netzer 



20061 ). If this is 



of random motions of emitters (as e.g. in extended envelopes, see 
the case, it might be that Fe II lines are also formed in the gas located in (or around) the 
starburst regions, i.e., the intensity of Fe II lines does not depend on the geometry of the 
BLR. The star forming regio ns can affect observed spectra of AGNs. As an example, recall 



here the results obtained by 



CroometaL 



( 120021 ). who found a strong anti-correlation 



with luminosity for the equivalent widths of [O II]A3727 and who suggested that the [O II] 
line observed in the high- luminosity AGNs may be emitted, to a large part, by intense 
star- forming regions; i.e., the AGN contribution to this line could be weaker than previously 
assumed. 

On the other hand, analysis of the AGN and AGN1 samples shows that there is no 
correlation between the widths (kinematics) of broad FL3 and luminosity. This would be 
expected if the widths were caused by rotational component of motion (by geometry), as is 
assumed in the classical BLR. Additionally, the EW Fe II shows a trend of being weaker as 
widths increase. This also can be caused by geometry, or it could be that as we look deeper 
into the BLR, the Fe II emission becomes weaker. 

In the case of a pure (or dominant) AGN emission there is a point-like photoionization 
source that influences the NLR emission: as the central source gets stronger, the continuum 
also becomes stronger, and it affects the processes in the NLR, as well as the size of the 
region. The Fe II emission becomes stronger. This may be caused by additional atomic 
processes in some parts of the BLR. In this case, the BLR is formed around the central 
black hole, and its geometry and gravitational motion influences line profiles (widths). 
Since dimensions will be affected by the central source luminosity there will be correlations 
between the broad line EWs and luminosity. However, the BLR geometry can be quite 
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different and one cannot expect correlations between the broad line FWHM (FWM10%) 
and luminosity. It is clear that in this case, the Baldwin effect and BG92 EV 1 are caused 
by the luminosity of central so urce and line-forming p rocesses in the NLR and part of the 



BLR (probably in the ILR, see 



Kovacevic et al. 



2010). 



3.1.2. Accretion rate 



The ratio between the bolometric and Eddington luminosities can be taken as equal 
to the dimensionless accretion rate m ~ L/M and it is well known that the Boroson and 



Green EV1 is closely related to the Eddington ratio, L/M (see e.g. 



Marziani et al. 



Boroson 



1985 



200 2]). The mass of a black hole (M) can be estimated using ( Wandel fc Yahil 



Padovani fc Rafanelli 



1988 



Vestergaard & Peterson 



2006 



McGill 



2001 



20081 ) 



M ~ Llf 00 (FWHM) 2 

where FWHM is the full width at half maximum of the broad component. From this, one 
can extract an expected relationship between the luminosity, accretion rate and FWHM as 



log(L) = const + 2 • log(L/M) + 4 • \og(FWHM). 

Consequently, if two groups of AGNs (in our case SB and SB1) have similar Eddington 
ratios, there will be a high correlation between the luminosity and FWHM. It may be an 
explanation for a good correlation between luminosity and FWHM seen in the SB and SB1, 
but the question remains, why in the AGN and AGN1 subsamples is there no correlation? 
Also, why is this correlation connected with the ratio of the narrow lines? One solution 
might be that SB and SB1 are young objects where starburst regions also contribute to the 
line emission (especially to the narrow lines), and that they have a similar accretion rate. 
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In that case we can expect a good correlation between the L and FWHM, and dominant 
starburst emission in narrow lines. 



Conclusions 



In this paper we have investigated connections between the Baldwin effect and the 
BG92 EV1 and different correlations between spectral characteristics of a sample of 302 
AGNs; dividing the sample with respect to the narrow [O III] and H/5 ratio and also by 
using the ratio of widths of [O III] and H/3. 

We have investigated different correlations within the subsamples. We found that: 



i) It seems that the line width of 
(FWHM=3000 km/s or 4000 km/s, see 



2010 



Kovacevic et al. 



2010 



jroad line Sv' 


Grupe 


2004 


. i 



w some authors 



Sulentic et al. 



2009 



Zamfir et al. 



etc.) is not relevant for separation of AGNs with broad lines 
into two groups (BLS1 and NLS1, see Fig. [6]), since some AGNs with FWHM< 3000 km/s 
or < 4000 km/s have the same characteristics as AGNs with broader lines. We propose 



to divide the subsamples of AGNs using as a criterion the flux ratio of the narrow [O III 
and narrow H/3 lines, or by using 1 
narrow and H/3 broad components 



and narrow H/3 lines, or by using the ratio of the FWHMs of [O III] and broad H/3, or H/3 



ii) There is a significant difference in the correlations between line properties for the 
two subsamples. This indicates that there is a significant difference in the physics of 
emitting gas and the origin of the broad line components. We found that BG92 EV1 and 
the Baldwin effect have the same physical background in both subsamples. In the SB and 



3 Since the width of the H/3 narrow and the corresponding [O III] line component is 
essentially the same 
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SB1 subsamples, eigenvector 1 shows strong acorrelation between luminosity, the FWHM 
of broad H/3, and the EWs of lines, while in the other subsamples there are high correlations 
only between luminosity and EWs of lines. 

hi) The narrow and broad line widths of the SB and SB1 subsamples are luminosity 
dependent and they seem not to be connected with only predominant rotational motion, 
but rather by randomly distributed high velocity gas. Conversely, in the case of the AGN 
and AGN1 subsamples, there is no FWHM^ broad line dependence on luminosity and the 
widths seem to be influenced by the BLR geometry. Such a high correlation between the 
luminosity and FWHM in the SB and SB1 subsamples may be explained as: a) a connection 
between starburst regions and the region where the broad H/3 line is formed, and b) the SB 
and SB1 subgroups having very similar accretion rates, in which case one can expect a high 
correlation between the luminosity and FWHM. 
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A. Estimation of the narrow H/3 component 

As we mentioned above, estimating the narrow H/3 component is a problem, especially 
in the case of broad-line AGNs. In principle, due to complex kinematics in the NLR the 
narrow lines usually cannot be represented by a pure Gaussian profile. But, if one would like 
to estimate the contribution of the NLR to the whole line, a Gaussian profile approximation 
may be used. This also gives some indication about velocities in the subregions as well 
as the contribution of the NLR to the total H/3 flux. Furthermore there is the problem 
of non- uniqueness in the solutions, i.e., that combinations of different Gaussians can give 
different contributions to the total H/3 line. 

In order to estimate the error-bars of the fits of the H/3 NLR intensity we perform the 
following procedure: we started from the best fit obtained, and we considered the different 
deviation of the H/3 NLR intensity as / x Ih/3nlr, where Ih/3nlr is the H/3 NLR intensity 
obtained from the best fit. The / values are varied from 0.9 to 0.1 and from 1.1 to 1.9 with 
a step of 0.1. For each value of / x Ih/3nlr we again looked for the best fit and found values 
of x 2 - We then constructed xVXo vs - R diagramsj (see Figs. [TOl and [TTT) . We found three 
characteristic X 2 /Xo vs - R shapes: (i) with a deep minimum, (ii) with a broad minimum, 
and also (iii) with an almost flat minimum in a wide interval of the narrow H/3 intensity 
(see Figs. [TOT and [TTI) . Next we assumed that a change of 5% of the \ 2 ma y be taken as an 
indicator of the error-bar, i.e., x 2 /Xo = 1-05 (horizontal lines in Figs. [TUl and [TT]) . 

In additional, we inspected spectra in the H/3+[OIII] spectral region, and found that 
there are three characteristic spectra (as shown in Fig. [TOl : (i) where the narrow H/3 is 
prominent and easy to fit, (ii) where the narrow H/3 is not so prominent, but can still 
be clearly seen and (iii) where the narrow H/3 cannot be clearly seen, but the best fit 



4 Xo is the value obtained in the primary fit 
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indicate existence of the narrow component. As can be seen in Fig. [10] the [O III] lines 
usually can be fit well with two Gaussians for each [O III] line, but the fitting procedure of 
the H/3 NLR line could sometimes be very difficult because the narrow component of H/3 
cannot be distinguished well especially in the third case. For that reason, the error in the 
determination of the H/3 NLR intensity dominates the error of the [O III]/H/3 NLR ratio. 
In Fig. [TUwe present three cases of the spectra where the narrow H/3 component is too 
weak, and, as can be seen, there are also cases where the y 2 test shows a very important 
contribution of the narrow H/3 component. 

Using this test we found 84 cases from the AGN and one case from the SB subsample, 
where the H/3 NLR intensities changed in a broad range, but x 2 shows only very small 
changes, which cannot reach 5% of Xo- This means that the contribution of the narrow 
H/3 in these objects is negligible. We also found 26 objects in the SB subsample for which 
the R =[0 III]/H/3 ratio has an error-bar which indicates that it could be higher than 
0.5. Finally, in the SB subsample we found 69 objects with R + AR clearly < 0.5. It is 
interesting that the correlations obtained are not changed when we exclude the spectra 
where the estimates of the narrow H/3 contribution have large error-bars. 

B. Error-bars of the FWHM and FWHM10% estimates 

As also mentioned above, the H/3 broad component was taken to be the sum of the H/3 
ILR and H/3 VBLR components (see Fig. [2]). Since there may be a problem in subtraction 
of the narrow component, we did a test by considering first the objects with flat minimum 
(i.e., where the narrow H/3 component cannot be properly estimated). We then fixed 
the narrow component taking / = 0.1 and / = 1.9 and found the best fit of the broad 
component. We measured FWHM and FWMI10% from these two broad components (using 
the ILR and VBLR) and found that the error-bars were in the range of 5%-7% (in all cases 
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less than 10%) in the case where the FWHM has a value between 1000 and 2000 km s 
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Fig. 1. — The best fit of the H/3 region (see text). 
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Fig. 2. — Measurement of the full width at half maximum and at 10% of the maximum (see 
text). 
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Fig. 3. — The BPT diagram for the subsample of 137 AGNs containing the narrow H/3 and 
Ha lines. Open circles denote AGNs with FWHM of the broad H/3 smaller than 3000 km 
s _1 , while full circles represent AGNs with FWHM > 3000 km s _1 . Arrows indicate upper 
limits for the AGN where maximal [N II] emission has been estimated (since the line was 
too weak). The horizontal dashed line represents the assumed border between 'AGN' and 
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Fig. 4. — The continuum luminosity vs. FWHM of the narrow [O III] lines. The left panel 
is for R < 0.5, and the right panel for R > 0.5 
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Fig. 6. — The continuum luminosity vs. the ratio of FWHMI [O III] and broad H/3. 
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FWM10% H/3 (lower panels) for AGNs with R <0.5 (left panels) and R >0.5 (right panels). 
Notation (see Appendix): full squares - R + AR <0.5, open squares - R + AR >0.5, crosses 
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Fig. 8. — The same as in Fig. [7] but for EW Fe II. 
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Fig. 9. — Continuum luminosity against FWHM of the [0 III] and broad H/3 lines for 
i?i > -0.8 (left) and R x < -0.8 (right). 
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Fig. 10. — Three characteristic spectra (left) and corresponding x 2 /Xo (right) against R. 
The first panel shows the spectra where narrow H/3 can be clearly seen. The second shows 
where it is present but covered by the ILR component; and the third shows the case where 
the narrow line cannot be clearly seen. The horizontal line denotes x 2 /Xo = l-05 and vertical 
line R=0.5. 
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Fig. 11. — The same as in Fig. [101 but for objects where the narrow H/3 is not prominent: 
(a) objects from the group with flat minima of x 2 {x 2 changes are less then 5% for different 
values of the H/3 NLR intensity), (b) objects with \ 2 changes higher than 5%, but less then 
10%, and (c) for objects with the very deep minimum of % 2 for different values of the H/3 
NLR intensity (more than 10%); 



Table 1: Spearman rank order correlations for the total sample and for subsamples: (1) log([0 III]5007/H/3 NLR)> 
and (2)log([0 III]5007/H/3 NLR)< 0.5 . 
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0.29 


2.3E-5 


0.88 





-0.05 


0.51 


-0.38 


1.4E-8 


-0.25 


2.4E-4 


0.42 


1.5E-10 




(i) 




0.81 





0.89 





-0.48 


1.5E-6 


0.27 


0.01 


-0.21 


0.05 


0.49 


6.4E-7 




(2) 



Table 2: Spearman rank order correlations for the total sample and for subsamples: (1) log(FWHM [0 
III]narrow/FWHM H/3)< -0.8 and (2)log(FWHM [O III]/FWHM H/3)> -0.8 . 





log(Ai 510 o) 


log(FWHM H/3) 


log(EW [O III]) 


log(EW Fc II) 


log(EW H/3 NLR) 


log(EW H/3 broad) 


log(FWMI10%H/3) 




r 


P 


r 


P 


r 


P 


r 


P 


r 


P 


r 


P 


r 


P 




Total sample 


0.42 


2.3E-14 


-0.45 


2.2E-16 


0.27 


1.6E-6 


-0.41 


9.4E-14 


0.14 


0.02 


0.43 


5.1E-15 


log (AL5100) 




(i) 


0.21 


0.01 


-0.45 


2.1E-8 


0.25 


0.003 


-0.56 


1.2E-12 


-0.08 


0.35 


0.23 


0.005 






(2) 


0.74 





-0.44 


2.5E-8 


0.32 


8.0E-5 


-0.37 


5.9E-6 


0.34 


3.1E-5 


0.74 







0.42 


2.3E-14 


total sample 


-0.07 


0.24 


-0.24 


2.3E-5 


-0.34 


1.5E-9 


0.48 





0.90 





log(FWHM H/3) 


0.21 


0.01 




(i) 


0.09 


0.29 


-0.41 


5.9E-7 


-0.12 


0.14 


0.32 


1.1E-4 


0.87 







0.74 







(2) 


-0.32 


1.0E-4 


0.24 


0.003 


-0.21 


0.009 


0.50 


1.6E-10 


0.87 







-0.45 


2.2E-16 


-0.07 


0.24 


total sample 


-0.41 


6.9E-14 


0.32 


2.1E-8 


0.24 


2.4E-5 


-0.05 


0.38 


log(EW [O III]) 


-0.45 


2.1E-8 


0.09 


0.29 




(i) 


-0.42 


3.0E-7 


0.71 





0.46 


1.7E-8 


0.13 


0.13 




-0.44 


2.5E-8 


-0.32 


1.0E-4 




(2) 


-0.43 


7.2E-8 


0.18 


0.03 


0.02 


0.81 


-0.33 


6.4E-5 




0.27 


1.6E-6 


-0.24 


2.3E-5 


-0.41 


6.9E-14 


total sample 


-0.01 


0.80 


0.07 


0.23 


-0.26 


6.3E-6 


log(EW Fc II) 


0.25 


0.003 


-0.41 


5.9E-7 


-0.42 


3.0E-7 




(i) 


-0.36 


1.1E-5 


0.05 


0.55 


-0.42 


2.0E-7 




0.32 


8.0E-5 


0.24 


0.003 


-0.43 


7.2E-8 




(2) 


-0.04 


0.65 


0.36 


6.7E-6 


0.18 


0.03 




-0.41 


9.4E-14 


-0.34 


1.5E-9 


0.32 


2.1E-8 


-0.01 


0.80 


total sample 


-0.16 


0.005 


-0.36 


1.6E-10 


log(EW H/3 NLR) 


-0.56 


1.2E-12 


-0.12 


0.14 


0.71 





-0.36 


1.1E-5 




(i) 


0.15 


0.08 


-0.15 


0.08 




-0.37 


5.9E-6 


-0.21 


0.009 


0.18 


0.03 


-0.04 


0.65 




(2) 


-0.19 


0.02 


-0.28 


6.8E-4 




0.14 


0.02 


0.48 





0.24 


2.4E-5 


0.07 


0.23 


-0.16 


0.005 


total sample 


0.50 





log(EW H/3 broad) 


-0.08 


0.35 


0.32 


1.1E-4 


0.46 


1.7E-8 


0.05 


0.55 


0.15 


0.08 




(i) 


0.38 


4.3E-6 




0.34 


3.1E-5 


0.50 


1.6E-10 


0.02 


0.81 


0.36 


6.7E-6 


-0.19 


0.02 




(2) 


0.46 


4.9E-9 




0.43 


5.1E-15 


0.90 





-0.05 


0.38 


-0.26 


6.3E-6 


-0.36 


1.6E-10 


0.50 





total sample 


log(FWMI10%H/3) 


0.23 


0.005 


0.87 





0.13 


0.13 


-0.42 


2.0E-7 


-0.15 


0.08 


0.38 


4.3E-6 




(i) 




0.74 





0.87 





-0.33 


6.4E-5 


0.18 


0.03 


-0.28 


6.8E-4 


0.46 


4.9E-9 




(2) 



-34 - 

Table 3: Principal component analysis of subsamples dependent of [O III]/H/3jvl.r and total 
sample. 



loa-ffO TTIl/H/^Arr o >0 5 subsamnle 


Comn 1 


Comn 2 

V. Will 1 7 . — 


Comn 3 


Standard deviation 


1 Kin 


l.zo f 


1.U40 


Proportion of Variance 


0.385 


0.255 


0.182 


Cumulative Pronortion 


0.385 


0.640 


0.822 


A-Zv5100 


-0.523 


0.088 


0.009 


FWHM H/3 


-0.154 


0.737 


-0.135 


ew [0 ml 


0.571 


0.138 


0.119 


EW Fe II 


-0.307 


-0.382 


0.676 


EW H/3 NLR 


0.529 


-0.124 


0.158 


EW H/3 broad 


0.057 


0.518 


0.696 


loaffO IIIl/H/3wr r? <0 5 subsamnle 

IWc \ V / llin J. l-i-J jY Li IX \J • *J OH U O tl 111 I J 1 1^ 


Comn 1 


Comn 2 


Comn 3 

v.. 1.1.1. 1/ . 


Standard deviation 


1. / Uo 


1 1 C7 


n nno 


Proportion of Variance 


0.485 


0.223 


0.136 


Cumulative Pronortion 


0.485 


0.708 


0.844 




-0.512 


-0.044 


-0.377 


FWHM H/3 


-0.490 


-0.111 


-0.446 


ew [o ml 

J i V T V ./ 1111 


0.384 


-0.532 


-0.160 


EW Fe II 


-0.331 


-0.326 


0.758 


EW H/3 NLR 


0.277 


-0.613 


-0.225 


EW H/3 broad 


-0.404 


-0.469 


0.091 


UU LOil OC1111JJ1C 


Comn 1 


Comn 2 


Comn 3 


Standard deviation 


1 /ICO 

1.4oo 


1.2 1 1 


1 noo 
l.Uzo 


Proportion of Variance 


0.354 


0.270 


0.174 


Cumulative Proportion 


0.354 


0.624 


0.798 


A-^5100 


-0.570 


-0.067 


0.004 


FWHM H/3 


-0.394 


0.544 


-0.117 


EW [0 III] 


0.454 


0.458 


0.178 


EW Fe II 


-0.246 


-0.471 


0.651 


EW H/3 NLR 


0.459 


-0.090 


0.363 


EW H/3 broad 


-0.208 


0.510 


0.631 
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Table 4: Principal component analysis of subsamples dependent of log(FWHM 
[O III]/FWHM H/3) ratio. 



lofffFWHM fO III1/FWHM H/9W -0 8 

1 V 1 _L.I_.LV_L 1 v__» _L_L_L / _L VV _L_L_LV_L _L_LA_y J \J .KJ 


m i h^p tyitYIp 

O U U O CI 1 1 1 IV 1 \j 


Comn 1 


Comn 2 

V_. Will 1 J . — 


Comn 3 

V_. v_/lll J V . -.7 


m~ p n H p tH H Pin prion 

O LjOj11U.cH U. LLC V IclllUll 




1 531 


1 .ZjUO 


1 07D 


Prnnnrtinn of V_P.ri_P.npp 

J. 1 UUUl U1U11 Ul V dl lallLL 




f) 3Q1 


0.241 


191 


Cumulative Proportion 




u.oyi 


n coo 


U.OZO 


^-^5100 




U.TlO 1 


-fl 948 




FWHM WR 

1 VV J-J-IVJ- ±±IJ 






-fl 7fi3 


-fl 194 


F,W fO TTTl 

1_J V V w 111 




n 585 




fl 1 38 


F,W Fp TT 

J__j VV -T t; 11 




u.ouz. 


99^ 


u.uou 


FW H/9 NFR 






fl 981 


116 


JtLiW rip Droaa 




n onn 


-U.4 < 1 


U.O ( 1 


log(FWHM [0 III] /FWHM H/3)> -0.8 


subsample 


Comp.l 


Comp.2 


Comp.3 


Standard deviation 




1.577 


1.047 


1.004 


Proportion of Variance 




0.414 


0.183 


0.168 


Cumulative Proportion 




0.414 


0.597 


0.765 






-0.530 


-0.014 


0.245 


FWHM H/3 




-0.530 


0.137 


0.065 


EW [0 III] 




0.302 


0.763 


-0.083 


EW Fe II 




-0.365 


-0.265 


-0.514 


EW H/3 NLR 




0.266 


-0.250 


-0.664 


EW H/3 broad 




-0.378 


0.516 


-0.472 



